Search for Lorentz and CPT Violation Effects in Muon Spin Precession 
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The spin precession frequency of muons stored in the {g — 2) storage ring has been analyzed for 
evidence of Lorentz and CPT violation. Two Lorentz and CPT violation signatures were searched 
for: a nonzero Auia {~^^a ~ )i and a sidereal variation of tj^ . No significant effect is found, 
and the following limits on the standard-model extension parameters are obtained: bz ~ —(1.0 ± 
1.1) X 10"^^ GeV; (m^dzo + Hxy) = (1.8 ± 6.0 x 10~^^) GeV; and the 95% confidence level limits 

2" GeV and fo'f" < 2.6 x 10"^'' GeV. 



< 1.4 x 10 



PACS numbers: ll.30.Cp, ll.30.Er, 13.40.Em, 12.20.Fv, 14.60.Ef 



The minimal standard model of particle physics is 
Lorentz and CPT invariant. Since the standard model 
is expected to be the low-energy limit of a more fun- 
damental theory such as string theory that incorporates 
gravity, Lorentz and CPT invariance might be broken 
spontaneously in the underlying theory flj . At low ener- 
gies, the Lorentz and CPT violation signals are expected 
to be small but perhaps observable in precision experi- 
ments. 

To describe the effects of spontaneous breaking of 
Lorentz and CPT invariance, CoUaday and Kostelecky [2] 
proposed a general standard model extension that can 
be viewed as the low-energy limit of a Lorentz covariant 
theory. Lorentz and CPT violating terms are introduced 
into the Lagrangian as a way of modeling the effect of 
spontaneous symmetry breaking in the underlying fun- 
damental theory. Other conventional properties of quan- 
tum field theory such as gauge invariance, renormaliz- 
ability and energy conservation are maintained, and the 



effective theory can be quantized by the conventional ap- 
proach. In a subsequent paper, Bluhm, Kostelecky and 
Lane discussed specific precision experiments with muons 
that could be sensitive to the CPT and Lorentz violating 
interactions 

In this letter we present our analysis for CPT and 
Lorentz violating interactions in the anomalous spin pre- 
cession frequency, uja, of the muon moving in a mag- 
netic field. In experiment E821 Q at the Brookhaven 
National Laboratory Alternating Gradient Synchrotron, 
muons are stored in a magnetic storage ring that uses 
electrostatic quadupoles for vertical focusing. The stor- 
age ring has a highly uniform magnetic field with a 
central value of Bq — 1.45 T, and a central radius of 
p = 7.112 m. Polarized muons are injected into the stor- 
age ring, and the positrons (electrons) from the parity- 
violating decay /i"*"*^") — > e"*"^"^ ^m(^m) ^e(^e) carry aver- 
age information on the muon spin direction at the time 
of the decay. Twenty-four electromagnetic calorimeters 
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around the ring provide the arrival time and energy of 
the decay positrons. As the muon spin precesses rela- 
tive to the momentum with the frequency oja, which is 
the the difference between the spin precession frequency 
ujs and the momentum (cyclotron) frequency ojc [J| , the 
number of high-energy positrons is modulated by uja- In 
the approximation that /3 • B = 0, 
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where 7 = l/y^l — (3'^. The anomaly a^i is related to the 
spin g-factor by = (g^ — 2)/2, with the magnetic mo- 
ment given by /Z = g{q/2m)s^ and q = ±e. The 'magic' 
momentum of Pm = 3.09 GeV/c, (7™ = 29.3), was used 
in E821 so that the second term in Eq. [1] vanishes, and 
the electric field does not contribute to uia- 

The magnetic field is measured using nuclear magnetic 
resonance (NMR) techniques [1]. The NMR frequency, 
averaged over both the muon distribution and the data 
collection period, is tied through calibration to the Lar- 
mor frequency for a free proton, and is denoted by ojp Q . 
Thus two frequencies are measured, cua and ujp. 

In the analysis presented here, the muon frequency uja 
is obtained from a fit of the positron (electron) arrival- 
time spectrum N{t) to the 5-parameter function 



N{t) = A^'oe"^ {l + Acos{ujat + (j))) 



(2) 



The normalization TVq, asymmetry A and phase (f> de- 
pend on the chosen energy threshold E. While more 
complicated fitting functions are used in the analysis 
for the muon anomaly _4], they represent small devia- 
tions from Eq. [2] and are not necessary for the CPT 
analysis. The anomalous magnetic moment is cal- 
culated from = Ti/{X — TV), where TZ = lOa/^p and 
A = ^Ji^,/^Jip = 3.183 345 39(10) (±30 ppb) is the muon- 
proton magnetic moment ratio [5| . 

For the muon, the Lorentz and CPT violating terms in 
the Lagrangian are @ 
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where iD\ = id\ — qA\, and the small parameters a^, 
bfi, Hkx, CkX cLkX represent the Lorentz and CPT viola- 
tion. All terms violate Lorentz invariance with —a^ip^'^ip 
and — &k'0757'*'0 CPT odd; all the other terms are CPT 
even. In this model the conventional figure of merit 
'^g = \9fi+ - 9t^- \/9average IS zcro at leading order; how- 
ever effects on the anomalous spin precession frequency 
uja do exist in lowest-order [s^. The frequency uJa is pro- 
portional to the magnetic field and therefore to LOp, so 
the sidereal variation of 7^ = uja/ojp is analyzed, rather 
than Wa directly. 

To compare results from different experiments, it is 
convenient to work in the non-rotating standard celes- 
tial equatorial frame {X,Y,Z} The Z axis is along 



the earth's rotational north pole, with the X and Y axes 
lying in the plane of the earth's equator. The preces- 
sion of the earth's rotational axis can be ignored because 
its precession period is 26,000 years. In this frame the 
correction to the (standard model) muon anomalous pre- 
cession frequency in Eq. [T] is calculated to be 

cos x + 2(6^ cos ilt + b'^ sinrii)sinx; (4) 



bf =±^+ m^djo + IsjklHkl. (J=X, Y, Z) , (5) 



where x is the geographic colatitude {— 90°- latitude) of 
the experiment location. For E821, xbnl = 49.1°. The 
sidereal angular frequency is O = 2t: /Ts, where Tg « 
23 hours 56 minutes. Eq. |4] predicts two signatures of 
Lorentz and CPT violation: a difference between the time 

averages of uj^ and u)^ , and oscillations in the values 
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of 
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uj'-' 



at the sidereal angular frequency. 



The E821 muon (g — 2) data have been analyzed for 
these two signatures. Data from the 1999 run (/x"*") 0], 
2000 ;8] and 2001 (/x") were used. Since bounds 
on clock comparisons of ^^^Hg and ^'^^Cs d, [l3 place 
limits on the Lorentz-violating energy shifts in the pro- 
ton precession frequency (wp) of ~ lO"^'' GeV, any shifts 
in the NMR measurements are at the mHz level, which 
is negligible compared to the uncertainty in the ampli- 
tude of any sidereal variation in Wq. A feedback system 
based on the reading from NMR probes at fixed locations 
around the ring keeps the field constant to within 3 ppm. 
Both the proton frequency LUp and the muon frequency LUa 
are measured relative to a clock stabilized to a LORAN C 
10 MHz frequency standard[Tl|. via a radio signal. The 
LORAN C frequency standard is based on Cesium hy- 
perfine transitions with mp = 0, which are insensitive to 
orientation of the clock. The Michelson-Morley type ex- 
periment of Brillet and Hall fl^ establishes that the frac- 
tional frequency shift. A/// of the LORAN radio signals 
due to the earth's rotation is less than 10^^^, so at the 
level of the precision of this experiment, the measured 
value of LUp is independent of sidereal time. Moreover, 
since LORAN C is the frequency standard for oja as well, 
any sidereal variation in the LORAN C standard would 
cancel in TZ. 

Comparison of the time averages of coa over many 
sidereal days (Eq. [4|) gives Aua = {^^a ) ~ i^a ) — 
(46^/7) cos X- For measurements made at different x 
and/ or Wp, 
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The colatitudes for the E821 fX^ and ^ measurements 
were identical, and the slightly different values of Wp Q 



3 



N 22SI1.2 

^ 2291 

(M 2290.8 

3 2290.6 

S 2290.4 



N 6179.16 
6179.15 
K 6179.14 
^ 6179.13 
S 6179.12 
X 10 ■ 





. 1 . . . . 1 . . . . ' 1 , . . . 1 . . 


. , 1 . . 





1000 2000 3000 4000 5000 

xlO^ 

time (second) 


affntwatt 








1000 2000 3000 4000 5000 

3 

X 10 

time (second) 










)0 5000 

3 

X 10 

time (second) 



FIG. 1; Values of ujait), ujp{t) and 7^(^) = ^(t) from the 

2001 run. The uncertainty on each uja point is about 
20 ppm. 



can be neglected. The E821 result Q is A7e = -(3.6 ± 
3.7) X 10~^, which corresponds to 



-(1.0 ± 1.1) X 10"^^ GeV or 



.7 ±8.9) X 10 



-24 



(7) 



a factor of 22 improvement over the limit that can be 
obtained from the CERN muon {g — 2) experiment p^ . 

The second signature of Lorcntz and CPT violation 
would be a variation of TZ{t) with a period of one sidereal 
day. The data are collected in 'runs' of approximately 60- 
minute duration, each with a time stamp at the beginning 
and end of run. Data from each of these different time 
intervals are fitted to the 5-parameter function (Eq. [J) to 
obtain uja, with the center of the time interval assigned as 
the time of that interval. The average magnetic field, cOp, 
for that time interval is used to determine TZ{t). Fig. [1] 
shows iiJait), i-^pit) and TZ(t) as a function of time for the 
2001 data collection period. 

A sidereal variation of uJa can be written as 
Wa(wp(tj),ii) = KojpiU) + A{iCOs{Qti + (f)). Dividing by 



TZ{U) = K 



A, 



LOp{ti) 



cos{rtti 



(8) 



where K — Aa^/(l + a^) is a constant, and Aq is the am- 
plitude of the sidereal variation with the sidereal period 
2n/fl. Two analysis techniques were used to search for an 



oscillation at the sidereal frequency: a multi-parameter 
fit to Eq. [51 and the Lomb-Scargle test [ll|, a spectral 
analysis technique developed for unevenly sampled data 
such as those displayed in Fig. [TJ With evenly sampled 
data it reduces to the usual Fourier analysis. For the 
time series {hi} with i = 1, . . . A^, the normalized Lomb 
power at frequency lu is defined as 



X;j^iC0s2[w(tj -r) 



Y^'i=i{hi ~ h) s\n[Lu{ti ~ r)] 



(9) 



Ei=iSin^[^(i» -^)] 
where h, a and r are defined as: 
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In searching for a periodic signal, the Lomb power is cal- 
culated over a set of frequencies. For a single frequency, 
with no corresponding periodic signal, the Lomb power 
is distributed exponentially with unit mean. If M in- 
dependent frequencies are scanned, the probability that 
none of them are characterized by a Lomb power greater 
than z is (1 — e^^)*^, assuming there is no signal present. 
The significance (confidence) level of any peak in Pn{^) 
is 1 — (1 — e~^)*^, which is the probability of the frequen- 
cies being scanned giving a Lomb power greater than z 
due to a statistical fluctuation. A small value of this 
probability therefore indicates the presence of a signif- 
icant periodic signal. In the case of equally separated 
points, the number of independent frequencies is almost 
equal to the number of time values. The lowest indepen- 
dent frequency, fa, is the inverse of the data's time span 
and the highest is roughly {N/2)fo, but, because of the 
uneven time sampling, may be somewhat greater. More 
generally, the number of independent frequencies, which 
depends on the number and spacing of the points, as well 
as the number of frequencies scanned, can be determined 
by Monte Carlo simulation, using Eq. [Hlto fit for M. 

The frequency spectrum for the 2001 data is shown in 
Fig. [21 The Lomb power at the sidereal frequency is 3.4. 
The probability that this is inconsistent with the null 
hypothesis is negligible. The Lomb power distribution of 
scanned frequencies (for M = 3144) shown in the lower 
half of Fig. [21 is consistent with an exponential distribu- 
tion, indicating there is no significant time- varying signal 
in the data. The Lomb power spectrum reaches only a 
modest maximum of 9.93 at fmax = 8.4 x 10^^ Hz, and 
is not present in the 1999 or 2000 data sets. If there were 
no oscillation signal at fmax, then 14% of the time the 
measured power is expected to be that large. Moreover, 
only a few of the Lomb power spectra taken from data 
subsets reveal a peak at that frequency. 



4 



Lomb-Scargle test on the simulated 2001 data 
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FIG. 2: The Lomb-Scargle test on the 2001 data (top plot). 
The horizontal lines show the confidence level associated with 
each Lomb power. At the sidereal frequency the Lomb power 
is 3.4, corresponding confidence level less than 10""^%. The 
Distribution of the Lomb power of the scanned frequencies 
(bottom plot) is consistent with an exponential distribution, 
indicating there is no time varying signal in the real data. 



FIG. 3: The Lomb-Scargle test on the simulated 2001 data 
with no signal. The shaded area shows the maximum Lomb 
powers greater than that of the real data. The locations of 
the maximum Lomb powers are randomly distributed (bottom 
plot). 



The significance of a possible sidereal signal in the 
power spectrum was carefully studied with a Monte Carlo 
simulation. First, a large number of artificial time spec- 
tra were generated. The distribution of times of the data 
points is chosen to be the same as for the real data. The 
value of each data point was distributed randomly with a 
central value equal to an average over all the actual data 
values, while having a standard deviation equal to that of 
each individual data point. The Lomb-Scargle test was 
then applied to 10000 simulated data groups. The distri- 
bution of maximum Lomb power and the corresponding 
frequencies are shown in Fig. [31 

Next, sidereal oscillation signals of different amplitudes 
Aq_ were introduced into the simulated data, and the 
resulting spectra were analyzed with both the multi- 
parameter fit and Lomb Scargle test. Table [T] lists the 
signal amplitude required so that 95% of the simulated 
data sets yield a larger or P(lo) than the real data. 
The two analysis methods give consistent results. 

Several potential systematic effects were studied. Since 
the sidereal period is very close to one solar day, the 
Lomb-Scargle test was applied to the data to check 
for false sidereal variations that might be produced by 
diurnal temperature changes. The upper limits on the 
amplitude of a sidereal variation in ojp were 0.04, 0.03 
and 0.08 ppm for the 1999, 2000 and 2001 data sets 



data set 


MPF amplitude 


L-S amplitude 




(ppm) 


(ppm) 


1999 /i+ 


5.5 


5.2 


2000 At+ 


2.2 


2.0 


1999/2000 


2.2 


2.0 


2001 iT 


4.2 


4.2 



TABLE I: The signal amplitude in ppm needed for 95% of 
the simulated data to have larger An, or -P(f2), than that of 
the real data. MPF means multi-parameter fit, L-S stands for 
Lomb-Scargle. 



respectively, significantly smaller than the limit on 
presented above. Additional studies were carried out on 
the 2001 data set. The data were folded back over a 
four sidereal day time period, (i.e. modulo four sidereal 
days), and then analyzed. To search for systematic ef- 
fects, other sub-window time periods, where no variation 
was expected, were also used, e.g. 24 hours (solar day), 
or an arbitrary number of minutes. 

No significant sidereal variation in 7^, and hence in to a-, 
is found in the E821 muon [g — 2) data. The limits on 
Ail from the MPF in Table U give at the 95% confidence 
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level: 

b'f = y^(fo^+)2 + (fe^+)2 < 1.4 X 10-24 GeV, ^^^^ 
K = + < 2.6 X 10-2^ GeV. 

For the dimensionless figure of merit obtained by dividing 
by ruf, H], we obtain r^^ eee 2 sin x^l^ /m^ < 2.0 x 10-^3 
and r^^ < 3.8 x lO^^s^ which can be compared with the 
ratio of the muon to Planck mass, m^/Mp = 8.7 x lO"^!. 
Using the E821 and CERN fTi| values of 7^ for ^+ and 
/i^ along with Eq. [51 we find 

(m^f^^o + Hxy) = (1-8 ± 6.0 x lO^^^) GeV. (12) 

An experiment that searched for sidereal variation 
in transitions between muonium hyperfine energy lev- 
els [l^, obtained < 5.0 x lO^^^. Penning trap ex- 
periments with a single trapped electron obtained r'^ < 
1.6 x 10-21 [l^. 



We find no significant Lorentz and CPT violation sig- 
natures in the E821 muon (g — 2) data, which we in- 
terpret in the framework of the standard model exten- 
sion 0. Limits on the parameters are of the order 10"23 
to 10"24 GeV, with the dimensionless figures of merit 
^ 10-2^. These results represent the best test of this 
model for leptons. Both r'^^^ and r^^^ are much less 
than ra^/Mp^ so E821 probes Lorentz and CPT viola- 
tion signatures beyond the Planck scale. 
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